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The spatial and temporal relationships between stellar age, kinematics, and chemistry are a funda¬ 
mental tool for uncovering the physics driving galaxy formation and evolution. Observationally, 
these trends are derived using carefully selected samples isolated via the application of appropri¬ 
ate magnitude, colour, and gravity selection functions of individual stars; conversely, the analysis 
of chemodynamical simulations of galaxies has traditionally been restricted to the age, metallic - 
ity, and kinematics of ‘composite’ stellar particles comprised of open cluster-mass simple stellar 
populations. As we enter the Gaia era, it is crucial that this approach changes, with simula¬ 
tions confronting data in a manner which better mimics the methodology employed by observers. 
Here, we use the SynCMD synthetic stellar populations tool to analyse the metallicity distri¬ 
bution function of a Milky Way-like simulated galaxy, employing an apparent magnitude plus 
gravity selection function similar to that employed by the RAdial Velocity Experiment (RAVE); 
we compare such an observationally-motivated approach with that traditionally adopted - i.e., 
spatial cuts alone - in order to illustrate the point that how one analyses a simulation can be, in 
some cases, just as important as the underlying sub-grid physics employed. 
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1. Introduction 

An inherent assumption underpinning the analysis of galaxy simulations is that they provide 
an analogous framework to that of galaxies in nature, in the sense that the chemo-kinetics of the 
massive stellar-like particles in simulations can be compared directly with the individual stars in 
observational samples. Considering that the typical mass of the ‘composite’ stellar particles in cos¬ 
mological, hydrodynamical, simulations of Milky Way-like galaxies is of order 10 5±1 M 0 , rather 
than the ~10 ()±1 M 0 encountered in nature, it is not unreasonable to question the importance of this 
mismatch in ‘resolution’; graphically, this ‘mismatch’ can perhaps best be appreciated by referring 
to Fig [I]. Further, observers construct samples of stars based upon their respective and individual 
colour, apparent magnitude, chemistry, and gravity, while simulators are limited primarily to the 
integrated luminosities and chemistry of the coeval, open cluster-scale, simple stellar populations 
used to represent the ‘composite’ stellar particles. 
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Figure 1: Left panel: A real colour-magnitude diagram comprised of ~ 10 4 stars drawn from NGC 6397 [9]. 
Right panel: How the same ~ 10 4 stars are represented by ‘composite’ stellar particles within a simulation. 


In what follows, we present a pilot study for our long-term programme aimed at developing 
powerful, and flexible, tools for better transforming simulations into the observer’s plane. Here, 
we employ the synthetic colour-magnitude diagram package SynCMD [5] to replace composite 
stellar particles with individual stars in a self-consistent fashion, thereby allowing stellar samples 
to be drawn by an observer situated within a simulation using apparent magnitude and gravity cri¬ 
teria. The impact of this observationally-motivated selection function upon the inferred metallicity 
distribution of an analogous ‘solar neighbourhood’, compared with the traditional approach (i.e., 
selecting composite stellar particles spatially), will be assessed, and our future directions outlined. 
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2. Methodology 

The characteristics of the ~L* Milky Way-analogue analysed for this work (MaGICC-gl536) 
are detailed in [4]. The simulation was realised with the smoothed particle hydrodynamics code 
Gasoline [10], employing the MaGICC (Making Galaxies in a Cosmological Context) stel¬ 
lar feedback prescription [3]. To transform the simulation into the observer’s plane, we use the 
SynCMD (synthetic colour-magnitude diagram tool) developed by [5]. Using the age, metallicity, 
and initial mass function underpinning stars in simulations, SynCMD decomposes each particle 
into its individual constituent (or ‘synthetic’) stars (of order 10 5 ‘synthetic’ stars per ‘compos¬ 
ite’ stellar particle, corresponding to ~10 11 individual stars for MaGICC-gl536), by mapping the 
composite particles onto the Padova isochrones ( [1]; [2]). 

3. Results 

For this pilot study, we situate the observer at an analogous ‘solar neighbourhood’, much as 
we did in [4]. From that vantage point, we ‘view’ the simulated galaxy, first isolating composite 
stellar particles within d=10 kpc of the observer and which avoid the galactic plane (i.e., enforcing 
a galactic latitude restriction of |Z?| >20 ° from the observer’s position). Such galactic plane restric¬ 
tions mimic the typical constraint imposed upon ground-based, optically-selected, all-sky surveys, 
such as RAVE [8]. 

For each composite stellar particle passing the spatial cut, we employ SynCMD to distribute 
the underlying ~ 10 s individual ‘synthetic’ stars, in the appropriate numbers (weighted by the 
initial mass function), into each ‘bin’ of the associated colour-magnitude diagram (CMD). After 
taking into account the distance modulus of each such star, the relevant colour-apparent magnitude 
diagram, from the observer’s vantage point, results. The upper panel of Fig [2] shows the resulting 
(I,V—I) CMD for the aforementioned stars passing the initial spatial cut; colour-coding here is by 
density of stars populating each CMD bin, as noted in the legend. 

Having transformed our composite stellar particles into individual synthetic stars, each with 
colours, luminosities, and gravities drawn from the Padova isochrones, we are now in a position 
to apply observationally-motivated selection criteria to the simulation. For this pilot study, we em¬ 
ploy criteria not dissimilar to that of RAVE [8]. We will show the results of three experiments: (i) 
including all stars with an I-band magnitude in the range 9</<12, regardless of gravity, which is 
reflected in the white horizontal lines in the upper panel of Fig ||; (ii) including both the aforemen¬ 
tioned I-band selection, and a surface gravity cut aimed at isolating main sequence and sub-giant 
branch (MS+SG) stars (logg>3.5); and (iii) including both the aforementioned I-band selection, 
and a surface gravity cut aimed at isolating giant branch (GB) stars (logg<3.5). The bottom left 
and right panels of Fig [| correspond to the colour-apparent magnitude diagrams associated with 
(ii) and (iii), respectively. 

To assess the impact of our observationally-motivated selection criteria in a quantitative sense, 
for this pilot study we examine the resulting metallicity distribution function (MDF) for each of 
the aforementioned SynCMD (i.e., ‘synthetic’) experiments (referred to as (i), (ii), and (iii) in the 
previous paragraph), and contrast those with the MDF inferred from using only the ‘composite’ 
stellar particles. 
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Figure 2: Upper: CMD for the synthetic stars passing the spatial cut, colour-coded by density. The two 
lines represent our I-band selection cut. Bottom left: As for the upper, but including gravity cut for main 
sequence + sub-giants. Bottom right: As for the upper, but including gravity cut for giant branch stars. 


The resulting four MDFs are shown in Fig [3j the thick solid line corresponds to the sub-set of 
‘composite’ stellar particles passing the initial spatial cut, while the thin solid line corresponds to 
the case in which also an I-band cut has been applied. One can immediately appreciate that the dis¬ 
tributions, while similar globally, differ in detail; formally, the MDF inferred from the composite 
simulation particles is significantly less skewed, with significantly lower kurtosis, than the counter¬ 
part derived employing a RAVE-like I-band selection cut (skewness=—1.2 vs —1.5; kurtosis=1.4 vs 
2.5). The differences become even more obvious when further constraining the observational cuts 
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beyond just apparent magnitude, to isolate main sequence plus sub-giant branch (MS+SG) stars 
(blue-dashed histogram in Fig |3j) or giant branch (GB) stars (red-dashed histogram in Fig |3j). In 
those cases, the skewness ranges from —2.2 (MS+SG) to —1.3 (GB), while the kurtosis ranges from 
7.1 (MS+SG) to 1.7 (GB). Differences in the mean and variance also range from 0.1 to 0.2 dex, 
across the four experiments. 



Figure 3: Comparative metallicity distribution function of the composite particles from our simulated galaxy 
(thick black line) and results from the decomposition into individual stars made with SynCMD (thin black 
line). We have further decomposed this synthetic population into main sequence + sub-giants (MS+SG, blue 
dashed line) and giant branch stars (GB, red dashed line), according to the log# restrictions described in § 3. 


In and of themselves, these quoted moments of the MDFs might not be enlightening, but a 
comparison with the corresponding moments from the broad suite of simulations described by [6] 
proves interesting. In [6], five variants of the same simulated disk galaxy were presented, employ¬ 
ing different initial mass functions, different treatments of radiation pressure from massive stars, 
different treatments of metal diffusion, and different treatments of the thermal cooling timescale. 
These five simulations were identical in every other consideration; as Table 2 in [6] shows, the 
resulting ranges in ‘solar neighbourhood’ MDF skewness and kurtosis vary from —0.9 to —1.8 
(skewness), and 0.9 to 3.8 (kurtosis). In other words, the impact of simply how one observes a 
simulation (i.e., observationally-motivated, as we have done here, or simple spatial cuts, as has 
been done previously by essentially the entire field) can be as quantitatively important as any of 
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the sub-grid physics treatments within the simulations themselves! 

In some sense, the results shown here should not be shocking, and are already appreciated 
by the observational community. For example, for the bright apparent magnitude cut employed 
here, the main sequence and sub-giant branch stars all lie within 1—2 kpc of the simulation ob¬ 
server, while the intrinsically brighter giants probe much greater distances. Without an a posteriori 
matching volume-limited selection criterion for the two samples (MS+SG vs GB), the different ef¬ 
fective distances of the samples mean that different metallicities and ages are encountered due to the 
interplay between the radial+vertical metallicity gradients and the spatially varying age-metallicity 
relations. 

4. Summary 

We have transformed a cosmological hydrodynamical simulation of an ~L* Milky Way-like 
galaxy from the simulator’s space to the observer’s plane, making use of the synthetic colour- 
magnitude diagram tool. Having replaced ~10 6 ‘composite’ stellar particles with r^j 10 11 ‘syn¬ 
thetic’ stars, we have employed various spatial, apparent magnitude, and surface gravity cuts to 
the individual synthetic stars, and viewed the simulation much as an observer would do in nature. 
This pilot study was not meant to be exhaustive, but instead, to illustrate that under very common 
conditions, selecting stars based upon actual observable characteristics, rather than just using mas¬ 
sive composite simulated stellar particles, can lead to inferred metallicity distributions which differ 
by as much as any of the various sub-grid physics treatments employed by simulators. While an 
enormous literature and community has evolved which is predicated upon understanding the sub¬ 
tleties of these sub-grid physics implementations, we have shown that it is equally as important 
to simply view the simulation ‘correctly’ as it is to get the sub-grid physics ‘right’. In the next 
phase of our work, we will also include the effect of foreground reddening, as well as explore 
the impact of magnitude-limited vs volume-limited sample selection upon not only the MDF, but 
also various chemistry-chemistry projections, the age-metallicity relation, and metallicity gradient 
determinations. 
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